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Abstract 

The most general exclusion single species one dimensional reaction-diffusion 
models with nearest-neighbor interactions which are both autonomous and 
can be solved exactly through full interval method are introduced. Using 
a generating function method, the general solution for, f„, the probability 
that n consecutive sites be full, is obtained. Some other correlation func- 
tions of number operators at nonadjacent sites are also explicitly obtained. 
It is shown that for a special choice of initial conditions some correlation 
functions of number operators called full intervals remain uncorrelated. 
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1 Introduction 



Most of the analytical studies on non-equilibrium statistical models belong to 
the low-dimensional (specially one dimensional) models fTkl3|. Analyzing one 
dimensional models, which are usually easier to investigate, helps us gaining 
knowledge on systems far from equilibrium. One of the techniques used to ob- 
tain exact results is the empty interval method (EIM), or its equivalent the full 
interval method (FIM). Among other things, it has been used to analyze the one 
dimensional dynamics of diffusion-limited coalescence [T4Hl7j . In these, one di- 
mensional diffusion-limited processes have been studied using EIM. There, some 
of the reaction rates have been taken infinite, and the models have been worked 
out on continuum. For the cases of finite reaction-rates, some approximate 
solutions have been obtained. Using this method £"„ (the probability that n 
consecutive sites be empty) has been calculated. (Alternatively, Fn is the prob- 
ability that n consecutive sites be full.) This method has been used to study 
a reaction-diffusion process with three-site interactions . EIM has been also 
generalized to study the kinetics of the q-state one dimensional Potts model 
in the zero-temperature limit [191 . In [20], all the one dimensional reaction- 
diffusion models with nearest neighbor interactions which can be exactly solved 
by EIM have been found and studied. Conditions have been obtained for the 
systems with finite reaction rates to be solvable via EIM, and then the equations 
of EIM have been solved. There solvability means that the evolution equation 
for En be closed. It turned out there, that certain relations between the re- 
action rates are needed, so that the system is solvable via EIM. When these 
conditions between reaction rates are met, the time derivative of E^s would be 
linear combinations of E^s. It was shown that if certain reactions are absent, 
namely reactions that produce particles in two adjacent empty sites, the coeffi- 
cients of the empty intervals in the evolution equation of the empty intervals are 
n- independent, so that the evolution equation can be easily solved. The criteria 
for solvability, and the solution of the empty-interval equation were generalized 
to cases of multi-species systems and multi-site interactions in f2TM23] . In [24] , 
models were studied which were solvable through EIM, but did include interac- 
tion which produce particles in two adjacent empty sites. There these models 
were investigated in continuum, although some terms in the evolution equation 
were missed, as will be discussed in the present paper. In [25], conventional EIM 
has been extended to a more generalized form. Using this extended version, a 
model has been studied which can not be solved by conventional EIM. Recently, 
the coagulation-diffusion process on a one dimensional chain has been studied 
using the empty- interval method [26| . There the behavior of the time-dependent 
double-empty- interval probability has been studied. In |27| . the exact two-time 
correlation and response functions for a one dimensional coagulation-diffusion 
process has been studied using EIM. 

In [28| , a ten-parameter family of reaction-diffusion processes was introduced 
for which the evolution equation of rt— point functions contains only n— or less- 
point functions. We call such systems autonomous. The expectation value of 
the particle-number in each site has been obtained exactly for these models. 
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In order to be an autonomous model there should be some constraints on the 
reaction rates. 

Here we study the most general exclusion single species one dimensional 
reaction-diffusion model with nearest-neighbor interactions, which can be solved 
exactly through the full interval method, and is autonomous. The reaction rates 
corresponding to these two models apart from corrections to [53], here lattice 
models are studied, while in |24j such models on continuum were investigated. 
The change from the empty interval to the full interval is, of course, not impor- 
tant; as a simple interchange of particles and holes would do that. The scheme 
of the paper is as follows. In Section 2, the most general exclusion single species 
one dimensional reaction-diffusion models with nearest-neighbor interactions are 
introduced, which are autonomous and can be solved exactly through FIM. In 
Section 3, the evolution equation for the full interval probabilities, F^s are 
obtained. In section 4, the steady state solutions for these probabilities are ob- 
tained, and then using a generating function method, the general solution for 
FnS is calculated. Correlation functions of number operators at nonadjacent 
sites are obtained in section 5. Finally, in section 6 some correlation functions 
of number operators are explicitly calculated for some special choice of initial 
conditions. These are probabilities of some disjoint parts of the lattice be full. 

2 Full interval and autonomy 

Consider a one dimensional lattice, any site of which is either occupied by a 
single particle or empty, and assume that the reactions, as well as the state of 
the system, are translationally invariant. Implicit in this, is that the lattice has 
no boundaries. But the lattice can still be finite, if it is circular. Defining Fn 
as the probability that n consecutive sites be full 



where an empty (occupied) site is denoted by by o (•), it is found (for example 
similar to [24]), that the most general single species nearest-neighbor interac- 
tions for which the evolution equations governing F„ 's are closed are 
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where ri,r2, qi, q2,w,wi and W2 are reaction rates. Then the Hamihonian for 
a two site interaction for models solvable through FIM is 
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The autonomy criteria leads to two more constraints ( [28] . for example) 

''1 + qi — Wi + w, 

r2 + q2 = W2 + w. (4) 
So an autonomous model solvable through FIM has five free parameters. 



3 The full interval evolution 

The full interval equation is 

dF. 

-T^ = {ri+ r2) [Fn-i + F„+i - 2 F„) - (qi + (72) {F,, - 
at 

- {n - I) [wi + W2 + w) Fn - [wi + W2 + 2 w) Fn+1, (5) 

where 

Fo := 1. (6) 

The difference of this with what obtained in (24j . apart from the obvious inter- 
change of particles and vacancies, is that the last term had been missed in [24] . 
Using (|4]), it is seen that the coefficient of Fn+i in the right-hand side of ([5]) 
vanishes. So, 

dF 

-^ = (7-1+ r-2) F„_i - [2 (n + r2) + qi + q2 + {n - 1) {wi + W2 + w)] F„. (7) 

Rescaling the time by {wi + W2 + w): 

i:= {wi+W2+w)t, (8) 
the evolution equation becomes 



dFn 
di 



where 



6F„_i - (a + n- 1)F„, (9) 
2 (ri -I- r2) -1-91-1-92 



W -I- Wl -I- W2 

n + r2 



111 + 101+102 



(10) 
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Equation Q, and the fact that the rates are nonegative guarantee that {wi + 
W2 + w) is positive, unless all of the rates are zero. These also show that 

a>l + b, (11) 
6 > 0. (12) 

From now on, the rescaled time i is denoted by t, so that the evolution equation 
is written as 

K=bF,,_i^{a + n-l)F„, (13) 
where dot means differentiation with respect to the rescaled time. 

4 The general solution for the full interval 

Let's first consider the large time values for F„ or the steady state solution. 
Denoting the time independent equation for the full interval by F^^, one has 

K'= J , K-i, (14) 

a + n — 1 

which combined with ([6]) results in 

C ^ (15) 

L [a + n) 

The general solution to is of the form 

n 

Fn{t) = + ^ c,„„ exp[-(a + m-l)t], (16) 



where c„m's are constants. Putting (1161) in (1131) . one arrives at 

{n - m) c„„i = fec„_i„i, m <n. (17) 
This allows one to determine c„m's in terms of Cmm's, which are denoted by 

dm- 



[n — m) ! 
(n — my. 



(18) 



So, 



^ ^n — m 

Fn it) = f^,f + V 7 TT drn exp[- (fl + m - 1) t] , (19) 

m=l ^ ' 



showing that the (largest) relaxation time r is obtained from 

1 



T=-. (20) 

a 
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The constants dm are to be obtained from the initial conditions. One way is to 
define the generating functions G and G^^ through 



Using these and 



G:=^i^„(0)x", 

oo 



n=0 



G{x) = G^'ix) — ^ ^ ' 



^-^ ^-^ (n — m)\ 

n=0 m=l ^ ' 



= G-(x) + ^ 5: 



' ^-^ in — m)! 

1 n—m 



d t'^ 



G^*(x) + ^d™ f-j exp(fea;), 

oo 

G"*(x) +exp(6a;) ^ d^x"". 



(21) 



(22) 



So, 



which resuhs in 



J2 d,nx"' ^ cxp(-fex) [G{x) - G^'^x)], 



E 

fc=0 



i-b) 



k r 



^ — k 



Via) 



r(a + m — k) 



For example, the occupation probability of any site is 



Fiit) = - 
a 



Fi{0)-- 
a 



cxp(— a t). 



(23) 



(24) 



(25) 



5 Correlation functions of number operators at 
nonadjacent site 

Denoting the number operator in the site i by n^, it is seen that 

(hi) - a (nj). (26) 



5 



This is in fact the same as the evolution equation for Fi, as it should be. The 
solution to p6)) is 

{n,){t) = (n,)(0) exp(-ai) + - [1 - exp(-at)]. (27) 

a 

Defining the correlation Ci-^...i^. as 

C^o-ifc := {nio---ntJ, (28) 
where no two indices are adjacent, it is seen that 

k 

C,,...,,^bY,C,^...r^...,^-{k + l)aa„....,, (29) 

where i means that the index i has been omitted. 

A simple change of variable makes the above equations simpler. Defining 

(again for the case no two indices are adjacent) one arrives at 

4-^. =-(fc + l)aC,„...... (31) 

The so called connected correlations C"^ are defined inductively through 

Q =: C^, (32) 

V 

where the summation runs over all partitions V of (jq • • -ik)- Such a partition 
divides the indices {io ■ ■ ■ ik) into a parts, where the /3'th part is denoted by 
Vfj^io ■ ■ ■ ik)- A simple induction shows that 

4...,^=-(fc + l)aC<;...,,. (34) 

Another induction shows that for fc > 0, adding a constant to any of n^ 's does 
not change the connected correlations. The reason is that, denoting the corre- 
lations and connected correlations corresponding to {ui^, + A) by a superscript 
A, one has 

where R contains terms in which the index j^. enters connected correlations of 
(to + 1) operators, where to is positive but less than k. Assuming that the 
independence of the connected correlations of (to + 1) operators is true for to 
positive and less than fc, it is seen that R does not depend on A. The sum of 
the first two terms of the right hand side is also obviously independent of A. So 
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the left hand side is independent of A as well. To complete the induction, one 
should prove that the connected correlation of two number operators does not 
change upon adding a constant term to each of them. This is obvious since in 
that case, R in the right hand side of (1551) is zero. 
In particular, one arrives at 

Cl....=Cl..,^, fc>0, (36) 

showing that 

C^.....W=C^o...,.(0)exp[-(fc+l)at], k>0. (37) 
The simplest example of this, is the connected two point function: 

{{nirij) - {ui) {nj)){t) = {{nirij) - {ui) {nj)){0) exp(-2at), > 1. (38) 

To consider correlations with possibly adjacent sites, let us first obtain the 
evolution equation for the full interval, without the assumption of translational 
invariance. Defining Fij as the probability that the sites beginning from i ending 
in j are full, one arrives at 



where 



F^j = hi F.j^i + 62 F,+ij -{a+j-i) F,j, (39) 



bi := 



W + Wi + W2 



b2 ■■= . (40) 

W + Wi + W2 

Then, reassuming translational invariance define Di as 

^= (n ( ft ne„+...+i,,_,+rn] ) , (41) 

\j=0 \m=l J I 

where 

£:=(4,...,^2fe), (42) 

and £j's are positive integers. This is the probability that £9 consecutive sites 
beginning from 1 be full, €2 consecutive sites beginning from (£0 + ^1 + 1) be 
full, and so on. In the special case fc = 0, this correlation is the same as the full 
interval: 

Di^Ft. (43) 

In the special case that all £j's with even j are one, D is reduced to the (fc + 1)- 
point function with nonadjacent sites: 

-Di,£i =Cio...i^, (44) 
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where 

im+l — im = ^2m+l + 1- 

Using and (HI]), it is then seen that 



(45) 



= Xl[^2£'£-e2,„+e2™-l + ^1 -D^-es ™+e2 ,„+i - (« + ''2 m - (46) 
m=0 

where is a (2 A: + l)-tuple, the only nonzero component of which is the m-th 
component being equal to one, and e_i and e2k+i are zero. 

A special case is the two-point function. For nonadjacent sites, one use 



Ci j — b (Ci + Cj ) 2 aCij, 



(47) 



the solution to which is 



Cij{t) = (n,,nj)(0) exp(~2 a i) + — [1 - cxp(-2 at)] 



(n,)(0) + (n,)(0) 



2b 



[exp(-a t) ~ exp(-2 a t)] , (48) 



where ([77]) has been used, and the fact that 

C,{t) - {n,){t). 

Using (p9)) . one has 

F,^- + -(a + l)i^,j, 

the solution to which is 

1.2 

-FjjXO = {n^ nj){Q) exp[-(a + 1) t] + 
+ 6i(n.)(0) + 62(n,)(0)- 



a (a + 1) 



{1 -cxp[-(a+ 1)<]} 



(49) 
(50) 



X {exp(— a t) — cxp[— (a + 1) t]}, j — i = 1- 
One thus arrives for an expression for the two point functions: 

rQ,(t), j-j>i 

(nj nj){t) = {F^ j{t), j - i^l . 

6 Uncorrelated full intervals 

A special case of the initial conditions is when 



(51) 



(52) 



^^(o)-n( n ^^o+-+^.,-i+™)(o), 

j=0 \m=l / 



(53) 
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which can be written as 

k 

Dm = X{F,,M^ (54) 

i=o 

Using it is seen that the ansatz 

k 

DE{t)^X{F,,^{t), (55) 

3=0 

does satisfy This does not mean that the system is completely uncorre- 

lated, as the ansatz 

F„(i) = (Fi )"(<), (56) 

does not satisfy But it means that the full intervals are uncorrelated to 

each other. Among other things, it does mean that the correlators C satisfy 

= {F^t+\ (57) 

if this holds initially. 
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